Gauge-Invariant Temperature Anisotropics and Primordial Non-Gaussianity 
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We provide the gauge-invariant expression for large-scale cosmic microwave background tempera- 
ture fluctuations at second-order in perturbation theory. It enables to unambiguously define the 
nonlinearity parameter /nl which is used by experimental collaborations to pin down the level of 
Non-Gaussianity in the temperature fluctuations. Furthermore, it contains a primordial term en- 
coding all the information about the Non-Gaussianity generated at primordial epochs and about the 
mechanism which gave rise to cosmological perturbations, thus neatly disentangling the primordial 
contribution to Non-Gaussianity from the one caused by the post-inflationary evolution. 
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Inflation has become the dominant paradigm to under- 
stand the initial conditions for the density perturbations 
in the early Universe which are the seeds for the Large- 
Scale Structure (LSS) and for the Cosmic Microwave 
Background (CMB) temperature anisotropies 0|. In the 
inflationary picture, primordial density and gravity-wave 
fluctuations are created from quantum fluctuations "red- 
shifted" out of the horizon during an early period of su- 
perluminal Universe expansion. Despite the simplicity of 
the inflationary paradigm, the mechanism by which cos- 
mological curvature (adiabatic) perturbations are gen- 
erated is not yet established. In the standard slow-roll 
inflationary scenario associated to one-single field, the 
inflaton, density perturbations are due to fluctuations of 
the inflaton itself when it slowly rolls down along its po- 
tential. In the curvaton mechanism Q the final curvature 
perturbation £ is produced from an initial isocurvature 
mode associated with the quantum fluctuations of a light 
scalar (other than the inflaton), the curvaton, whose en- 
ergy density is negligible during inflation. Recently, other 
mechanisms for the generation of cosmological perturba- 
tions have been proposed, the inhomogeneous reheating 
scenario Q , ghost-inflation and the D-cceleration sce- 
nario 5], to mention a few. A precise measurement of 
the spectral index uq of comoving curvature perturba- 
tions will provide a powerful constraint to slow-roll in- 
flation models and the standard scenario for the genera- 
tion of cosmological perturbations which predicts |nf — 1| 
significantly below unity. However, alternative mecha- 
nisms generically also predict a value of very close to 
unity. Thus, even a precise measurement of the spectral 
index will not allow us to efficiently discriminate among 
them. On the other hand, the lack gravity- wave signals in 
CMB anisotropies will not give us any information about 
the perturbation generation mechanism, since alternative 
mechanisms predict an amplitude of gravity waves far 
too small to be detectable by future experiments aimed 
at observing the B-mode of the CMB polarization. 

There is, however, a third observable which will prove 



fundamental in providing information about the mech- 
anism chosen by Nature to produce the structures we 
see today. It is the deviation from a Gaussian statistics, 
i.e., the presence of higher-order connected correlation 
functions of CMB anisotropies. Since for every scenario 
there exists a well defined prediction for the strength of 
Non-Gaussianity (NG) and its shape as a function of the 
parameters, testing the NG level of primordial fluctua- 
tions is one of the most powerful probes of inflation Q 
and is crucial to discriminate among different - but oth- 
erwise indistinguishable - mechanisms. For instance, the 
single-field slow-roll inflation model itself produces neg- 
ligible NG, and the dominant contribution comes from 
the evolution of the ubiquitous second-order perturba- 
tions after inflation, which is potentially detectable with 
future observations of CMB temperature and polariza- 
tion anisotropies. This effect must exist regardless of the 
inflationary models, setting the minimum NG level of 
cosmological perturbations. Therefore, if we do not find 
any evidence for this ubiquitous NG, then it will chal- 
lenge our understanding of the evolution of cosmological 
perturbations at a deeper level. 

Motivated by the extreme relevance of pursuing NG in 
the CMB anisotropies, in this Letter we provide the exact 
expression for large-scale CMB temperature fluctuations 
at second order in perturbation theory. This expression 
has various virtues. First, it is gauge-invariant. Second, 
from it one can unambiguously extract the exact defini- 
tion of the nonlinearity parameter /nl which is used by 
the experimental collaborations to pin down the level of 
NG in the temperature fluctuations. Third, it contains 
a "primordial" term encoding all the information about 
the NG generated in primordial epochs, namely during or 
immediately after inflation, and depends upon the var- 
ious fluctuation generation mechanisms. As such, the 
expression neatly disentangles the primordial contribu- 
tion to the NG from that arising after inflation. Finally, 
the expression applies to all scenarios for the generation 
of cosmological perturbations. 
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In order to obtain our gauge-independent formula for 
the temperature anisotropics we first perturb a spatially 
flat Robertson- Walker background. Here we follow the 
formalism of Ref. expanding metric perturbations in 
a first and a second-order part as 
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where the scale factor a(r\) is a function of the conformal 

time rj. The functions <j)^ r \ bj'f\ and Xij ■> where 
(r) = (1,2), stand for the rth-order perturbations of 
the metric. It is standard use to split the perturba- 
tions into the so-called scalar, vector and tensor parts, 
according to their transformation properties with respect 
to the 3-dimensional space with metric Sij , where scalar 
parts are related to a scalar potential, vector parts to 
transverse (divergence-free) vectors and tensor parts to 
transverse trace-free tensors. Thus <f> and ip are scalar 
perturbations, and for instance, uj^ — diUJ^ + u>l r \ 
where is the scalar part and uj^ is a transverse 
vector, i.e. crw 2 = 0. The metric perturbations will 
transform according to an infinitesimal change of coor- 
dinates. From now on we limit ourselves to a second- 
order time shift r\ — > rj — a^ + | (a'^a^ — &(2)), where 
a prime denotes differentiation w.r.t. conformal time. 
In general a gauge corresponds to a choice of coordi- 
nates defining a slicing of spacetime into hypersurfaces 
(at fixed time rj) and a threading into lines (correspond- 
ing to fixed spatial coordinates x), but in this Letter 
only the former is relevant so that gauge-invariant can 
be taken to mean independent of the slicing |2| • For ex- 
ample, under the time shift, the first-order spatial curva- 
ture perturbation transforms as ip^ — > tp^ —Ham 
(here H = a' /a), while ^« -> ^ + a' (1) + Ha^, 
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■ diaS x > , and the traceless part of the spatial 
metric y& turns out to be gauge- invariant. At second 
order in the perturbations we just give some useful ex- 
amples like the transformation of the energy density and 



the curvature perturbation 7] 5^ p — > S^ 2 
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In particular, there ex- 



The key point here is that the gauge-invariant comoving 
curvature perturbation C^ 2 ) remains constant on super- 
horizon scales after it has been generated and possible 
isocurvature perturbations are no longer present. There- 
fore, £( 2 ) provides all the necessary information about 
the "primordial" level of NG generated either during in- 
flation, as in the standard scenario, or immediately after 
it, as in the curvaton scenario. Different scenarios are 
characterized by different values of , while the post- 
inflationary nonlinear evolution due to gravity is common 
to all of them [EEIEIE]- For example, in standard 
single-field inflation is generated during inflation and 

its value is C (2) = 2 (C (1) ) 2 + O (n c - 1) USE!- 

We now construct in a gauge-invariant way tem- 
perature anisotropics at second order. Temperature 
anisotropics beyond the linear regime have been cal- 
culated in Refs. following the photons path from 
last-scattering to the observer in terms of perturbed 
geodesies. The linear temperature anisotropies read 
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where = + ( . . ._, v; ( 

script £ indicates that quantities are evaluated at last- 
scattering, e 1 is a spatial unit vector specifying the direc- 
tion of observation and the integral is evaluated along 
the line-of-sight parametrized by the affine parameter 
A. Eq. includes the intrinsic fractional tempera- 
ture fluctuation at emission T£, the Doppler effect due 
to emitter's velocity v^ 1 and the gravitational redshift 
of photons, including the Integrated Sachs- Wolfe (ISW) 
effect. We omitted monopoles due to the observer O 
{e.g. the gravitational potential ipQ evaluated at the 
event of observation), which, being independent of the 
angular coordinate, can be always recast into the defi- 
nition of temperature anisotropies [lfij ] . Notice however 
that the physical meaning of each contribution in Eq. @ 
is not gauge-invariant, as the different terms are gauge- 
dependent. However, it is easy to show that the whole 
expression is gauge- invariant. Since the temperature 
T is a scalar, the intrinsic temperature fluctuation trans- 



forms as 7g — > Tg i; + (T'/T)am = T y £ L ' — Ham, having 
used the fact that the temperature scales as T oc a -1 . 
Notice, instead, that the velocity v^- does not change. 
Therefore, using the transformations of metric perturba- 
tions we find 



ists an extension at second order of the well-known gauge- 
invariant variable C^ 1 -* = —tp^ — H^-^r- (the curvature 
perturbation on uniform density hypersurfaces). It is AT^ AT' 
given by C = C (1) + (V 2 )C (2) , where ~~ 
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where we have used the fact that the integral is evalu- 
ated along the line-of-sight which can be parametrized by 
the background geodesies x^ '^ = (A, (Xo — ^s)^) (with 



3 



d\/drj — 1), and the decomposition for the total deriva- 
tive along the path for a generic function /(A, x l (X)), 

f = Ia = lk + d ^ el - E( l' ® shows that the expres- 
sion for first-order temperature anisotropies is indeed 

gauge-invariant (up to monopole terms related to the ob- 
server O). Temperature anisotropies can be easily writ- 
ten in terms of particular combinations of perturbations 
which are manifestly gauge- invariant. For the gravita- 
tional potentials we consider the gauge-invariant defini- 
tions = and ^ = 0« + Hoj {1) + . 
For the (0 — i) component of the metric and the trace- f A d\ 



less part of the spatial metric we define oj. 
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For the matter variables we use a 



and x 

gauge-invariant intrinsic temperature fluctuation Tqj 
t' 1 ' — Huj^, while the velocity itself is gauge-invariant 
Vqj — v^ 1 under time shifts. Following the same 
steps leading to Eq. Q one gets the linear temperature 
anisotropies in Eq. (PJ in terms of these gauge-invariant 
quantities 
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where A™ = + Vgi ' ~< 2 Ay _ 

we omitted the subscript 5. For the primordial fluctua- 
tions we are interested in the large-scale modes set by the 
curvature perturbation £^ . Defining a gauge-invariant 
density perturbation 5^ pgi = 5^' p+p'cu^' , we write the 
curvature perturbation as £gi = — V'gi ~ H{5^ pGi/ p 1 )- 
Since for adiabatic perturbations in the radiation (7) and 
matter (m) eras (l/4)(5« Pl / Pl ) = (l/3)(<5« p m / p m ), 
one can write the intrinsic temperature fluctuation as 
T (D = (l/4)(<5( 1 V 7 / / 9 7 ) = -H{5^p/p') and a gauge- 
invariant definition is Tq{ — —Ti(S^pQi/p'). In the 
large-scale limit, from Einstein equations, in the mat- 
ter era 
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VGI — S^GI 



Thus we obtain the 

at' 1 ' 



large-scale limit of temperature anisotropies © 

2 V'gi + Cgi = V'gi/S, i-e. the usual Sachs- Wolfe effect. 

At second order, the procedure is similar to the one de- 
scribed so long, though more lengthy and cumbersome. 
We only provide the reader with the main steps to get 
the final expression. The second-order temperature fluc- 
tuations in terms of metric perturbations read [Tij 
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I\(X) are the photon wave vectors, with ii(A) 
given by the integral in Eq. and I{ (A) is obtained 
from the same integral replacing the time derivative 
with a spatial gradient. Finally in Eq. © x^°(X) = 
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and ^'(A) = 
-A')A«< 1 ] are 



the geodesies at first order, and — e 1 — is 
the direction of the photon emission. As usual we have 
omitted the monopole terms due to the observer. Us- 
ing the transformation rules of Ref. , it is possible to 
check that the expression @) is gauge- invariant. We can 
express the second-order anisotropies in terms of explic- 
itly gauge-invariant quantities, whose definition proceeds 
as for the linear case, by choosing the shifts such that 
yj{ r ) = 0. For example, we consider the gauge-invariant 
gravitational potential 01 
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where = ui^ + uj^oj^' + W~ 2 d^-A^ - 
2oj^ diUi( l >]. Expressing the second-order temperature 

anisotropies (jSJ in terms of our gauge-invariant quanti- 

(2) 

ties and taking the large-scale limit we find AT^'/T = 
(1/2)4^ - (1/2) (4V) 2 + (1/2)t$ + (having 

dropped the subscript £), and the gauge- invariant in- 

(2) 

trinsic temperature fluctuation at emission is r GI ' = 
(l/4)( ( j( 2 )p 7 J i/p 7 ) _ 3(4\ ) ) 2 . We have dropped those 
terms which represent integrated contributions and other 
second-order small-scale effects that can be distinguished 
from the large-scale part through their peculiar scale de- 
pendence. At this point we make use of Einstein's equa- 
tions. We take the expression for £( 2 ) in Eq. (J2J, and we 
use the (0 — 0) component and the traceless part of the 
(i—j) Einstein's equation at second order (see Eqs. (153) 
and (155) of Ref. Thus, on large scales we find that 
the temperature anisotropies are given by 
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kernel K, 



where we have defined 

iov- 4 a i ^'(5V (1) ^V' (1) ) - v~ 2 (f fltyw&v^)- 

Eq. (JSJ is the main result of this Letter. It clearly 
shows that there are two contributions to the final 
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nonlinearity in the large-scale temperature anisotropies. 
The contribution, [(^ - 2«£ I ) ) 2 ], comes from the 
"primordial" conditions set during or after inflation. 
They are encoded in the curvature perturbation £ which 
remains constant once it has been generated. The 
remaining part of Eq. (JSJ describes the post-inflation 
processing of the primordial non-Gaussian signal due 
to the nonlinear gravitational dynamics, including 
also second-order corrections at last scattering to the 
Sachs- Wolfe effect Thus, the expression in Eq. © 

allows to neatly disentangle the primordial contribution 
to NG from that coming from that arising after inflation. 
While the nonlinear evolution after inflation is the same 
in each scenario, the primordial content will depend on 
the particular mechanism generating the perturbations. 
We parametrize the primordial NG in the terms of 
the conserved curvature perturbation (in the radiation 
or matter dominated epochs) = 2a (C^) , where 
a depends on the physics of a given scenario. For 
example, in the curvaton case a = (3/4r) — r/2, where 
r w (At/p)d is the relative curvaton contribution to 
the total energy density at curvaton decay [§]. In 
the minimal picture for the inhomogeneous reheating 
scenario, a = 1/4. For the other scenarios we refer 
the reader to Ref.[f|. From Eq. JHJ) we can extract the 
nonlinearity parameter /nl which is usually adopted 
to phcnomenologically parametrize the NG level of 
cosmological perturbations and has become the standard 
quantity to be observationally constrained by CMB 
experiments 0, . The definition of /nl adopted in 
the analyses performed in Refs. 0,0] goes through the 
conventional Sachs- Wolfe formula AT/T — — <£>/3 where 
<f> is Bardeen's potential 0], which is conventionally 
expanded as (up to a constant offset, which only affects 
the temperature monopole) $ = $l + /nl * ($l) , with 
$l = —<Pgi- Here the * product reminds the fact that 
the nonlinearity parameter might have a non-trivial 



i) nam, 



scale dependence j6j. Therefore, using £W = — IV'gi 



during matter domination, from Eq. JSJ we read the 
nonlinearity parameter in momentum space 
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with 



where K = 10 (k x ■ k 3 )(k 2 • k 3 )/fc 4 - f k x • k 2 /fc 2 
k 3 + ki + k 2 = and k = |k 3 |. In fact the formula 
@ already accounts for an additional nonlinear effect 
entering in the CMB angular 3-point function from the 
angular averaging performed with a perturbed line ele- 
ment dfi(i - @, 

implying a +1 shift in /nl- In 
particular within the standard scenario where cosmolog- 
ical perturbations are due to the inflaton the primordial 



contribution to NG is given by a = 1 — | (n^ 
where the spectral index is expressed in terms of the usual 
slow- roll parameters as — 1 = — 6e + 2i] 1]. The non- 
linearity parameter from inflation now reads 



inf 
NL 



5 f 
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12 v c 
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Therefore the main NG contribution comes from the 
post-inflation evolution of the second-order perturbations 
which give rise to order-one coefficients, while the primor- 
dial contribution is proportional to — 1| <C 1. This 
is true even in the "squeezed" limit first discussed by 
Maldacena 19], where one of the wavenumbers is much 
smaller than the other two, e.g. k\ -C &2,3 an d /C — ► 0. 
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